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I. INTRODUCTION 

In theoretical cosmology, many insights can already be 
gained from spatially isotropic Friedmann-Robertson- 
Walker models 



ds 2 = -At 2 + a(r) 2 



1 — fcr 2 



+ r 2 (dtf 2 + sin 2 i?d<^ 2 ) 



. (1) 

with k = or ±1. The matter content in such a highly 
symmetric space-time can only be of the form of a perfect 
fluid with stress-energy tensor T a b = pu a Ub + P(g a b + 
u a Ub) where p is the energy density of the fluid, P its 
pressure and u a the 4- velocity vector field of isotropic co- 
moving observers. Once an equation of state P = P(p) 
is specified to characterize the matter ingredients, the 
continuity equation p + 3H(p + P) = with the Hubble 
parameter H = a/a allows one to determine the behavior 
of p(a) in which energy density changes during the expan- 
sion or contraction of the universe. This function, in turn, 
enters the Fricdmann equation H 2 + k/a 2 — 8nGp/3 and 
allows one to derive solutions for a(r). 

In general, one would expect the equation of state 
P = P{p) to be non-linear which would make an ex- 
plicit solution of the continuity and Friedmann equations 
difficult. It is thus quite fortunate that in many cases 
linear equations of state P — wp with w constant are 
sufficient to describe the main matter contributions en- 
countered in cosmology at least phenomenologically. The 
influence of compact objects on cosmological scales is, 
for instance, described well by the simple dust equation 
of state P(p) = 0. Relativistic matter, mainly electro- 
magnetic radiation, satisfies the linear equation of state 
P = ip. The latter example is an exact equation de- 
scribing the Maxwell field, rather than an approximation 
for large scale cosmology. It is thus, at first sight, rather 
surprising that the dynamics of electromagnetic waves in 
a universe can be summarized in such a simple equation 
of state irrespective of details of the field configuration. 
The result follows in the standard way from the trace- 



freedom of the electromagnetic stress-energy tensor and 
is thus related to the conformal symmetry of Maxwell's 
equations. That the availability of such a simple equation 
of state is very special for a matter field can be seen by 
taking the example of a scalar field <f> with potential V(<f>). 
In this case, we have an energy density p — ^<j) 2 + V(4>) 

and pressure P = \(\> 2 — V(4>). Unless the scalar is free 
and massless, V(<f>) — for which we have a stiff fluid 
P = p, there is no simple relation between pressure and 
energy density independently of a specific solution. 

Any conformal symmetry such as that of elecromag- 
netism might be broken by quantum effects especially 
when quantum gravity with its new scale provided by 
the Planck length is taken into account. The coupling of 
the electromagnetic field to geometry will then change, 
and exact conformal symmetries can easily be violated. 
Accordingly, one expects corrections from quantum grav- 
ity to the radiation equation of state and corresponding 
effects in the universe evolution during radiation domi- 
nated epochs. In loop quantum cosmology [1] equations 
of state of matter fields are in general modified by pertur- 
bative corrections at large scales and non-perturbative 
ones on small scales Q. This has mainly been stud- 
ied so far for a scalar field for which quantum modifi- 
cations can be so strong that negative pressure results 
independently of the chosen potential [3( . The main rea- 
son is the fact that the isotropic scalar field Hamiltonian 
H(f, = ^a~ 3 pi + a 3 V ((f)), where p<f, is the momentum of 
4>, contains an inverse power of the scale factor a. For 
quantum gravity, this factor has to be quantized, too. Us- 
ing the methods of 0], it turns out that inverse powers 
receive strong loop quantum corrections at small length 
scales Q • Accordingly, such modifications play a role for 
effective equations describing the universe after the big 
bang (or even during the quantum transition through the 
big bang singularity). During later stages, modifications 
are expected to decrease in size, but they might still be 
relevant due to sometimes tight constraints on evolution 
parameters. 
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An extension to the usual matter ingredients of cosmol- 
ogy with linear equations of state is, however, difficult 
since the modification is based on quantizations of the 
fundamental field Hamiltonians. Equations of state are 
obtained from fundamental Hamiltonians after an analy- 
sis of the matter field equations, which can be difficult in 
general especially when quantum effects are taken into 
account. The only exception is the dust case since it 
implies a constant Hamiltonian (the total mass of dust) 
which is straightforwardly quantized without any correc- 
tions. Thus, although the dust energy density is pro- 
portional to a~ 3 and metric dependent in a way which 
involves the inverse, it does not receive any modification 
since the Hamiltonian, i.e. total energy a 3 p, is the essen- 
tial object to be quantized. For radiation with p oc a~ 4 
the expectation is not clear since the total energy does 
behave like an inverse power of a, but this follows only af- 
ter an indirect analysis of the field dynamics. It is not the 
solution p(a) oc a~ 4 of the continuity equation which is 
quantized but the original field Hamiltonian from which 
the equation of state has to be derived first. One thus 
has to go back to the fundamental Maxwell Hamiltonian, 
derive energy density and pressure and see how quantum 
effects change the equation of state. If this is completed, 
one may attempt to solve the continuity equation to ob- 
tain corrections to p(a). 

We will derive such corrections in this article, using 
the canonical quantization given by loop quantum grav- 
ity [1, d, 0] • Candidates for Hamiltonian operators of the 
Maxwell field have been proposed [4( which show sev- 
eral sources of correction terms. To derive corrections to 
the equation of state, however, we need to perform the 
usual calculation in a Hamiltonian formulation. Thus, we 
first present the canonical formulation for the free classi- 
cal Maxwell field to rederive the standard result for the 
equation of state parameter w without reference to an 
action or the stress-energy tensor. Appropriate modifi- 
cations to the matter Hamiltonian Hm are then made to 
derive possible loop quantum gravity corrections to the 
equation of state w. We will show that one case of cor- 
rections results again in a linear equation of state, albeit 
in a corrected way which depends on the basic discrete- 
ness scale of quantum gravity. In this case we are able to 
express, as in the classical case, the full field dynamics in 
terms of a simple modified w, and to solve explicitly for 
p{a). Our derivation takes into account inhomogeneous 
field configurations and presents the first modified equa- 
tion of state obtained for a realistic matter source in loop 
quantum gravity. 

II. CANONICAL FORMULATION 

In a canonical formulation, the Hamiltonian Hm rather 
than the action is used to determine equations of motion 
of any function / on the phase space by means of Poisson 
brackets, / = {/, Hm}- The Poisson structure defines 



the kinematical arena which follows from the field vari- 
ables and momenta. The basic configuration variable in 
a Lagrangian formulation of Maxwell's field theory is the 
vector potential A a which determines the field strength 
tensor 

F ab = V a A b - V b A a , (2) 

where V a is the covariant derivative operator. Notice 
that V a can be replaced by the partial derivative op- 
erator d a even on a curved space-time since the field 
strength tensor F ab is antisymmetric. The action for the 
free Maxwell field in an arbitrary background space is 
given by 

S M = I d 4 x^F ab F ab 

16tt J 

= ~ f d 4 x^F ab F cd g ac g M (3) 
16tt J 

where g is the determinant of the Lorentzian space-time 
metric g ab - From the action one obtains Maxwell's equa- 
tions as the Euler-Lagrange equations extremizing Sm- 
A canonical formalism (Hamiltonian framework) is 
achieved by performing a Legendre transform of this ac- 
tion Sm, replacing time derivatives of configuration vari- 
ables by momenta. This, as always, requires one to treat 
space and time differently and is the reason why the 
canonical formulation is not manifestly covariant. We in- 
troduce a foliation of the space-time (M, g ab ) by a family 
of space- like hypersurfaces E t : t = constant in terms of a 
time function t on M. Canonical variables will depend 
on which time function one chooses, but the resulting 
dynamics of observable quantities will remain covariant. 
Furthermore, let t a be a timelike vector field whose inte- 
gral curves intersect each leaf E* of the foliation precisely 
once and which is normalized such that t a V ' a t = 1. This 
t a is the 'evolution vector field' along whose orbits differ- 
ent points on all E t = E can be identified. This allows us 
to write all space-time fields in terms of t-dependent com- 
ponents defined on a spatial manifold E. Lie derivatives 
of space-time fields along t a are identified with 'time- 
derivatives' of the spatial fields. 



A. Hamiltonian 

Let us, as illustrated in Fig. [TJ decompose t a into nor- 
mal and tangential parts with respect to E t by defin- 
ing the lapse function N and the shift vector N a as 
t a = Nn a + N a with N a n a = 0, where n a is the unit nor- 
mal vector field to the hypersurfaces E t . The space-time 
metric g ab induces a spatial metric q ab by the formula 
9ab = lab — n a n b . Now using n a = iV _1 (t a — N a ) and 

q ab _ g ab + n a n b 

to project fields normal and tangential 
to Et, we can decompose the field strength tensor F ab 
and the action Sm as follows: 
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FIG. 1: Decomposition of the evolution vector field t a in terms 
of the normal n a to spatial slices and a spacelike part N a . 



F ab n a = ^(F ab t a -N a F ab ) = ^(A b ~d b {A a t a )-N a F al 



F ab F ab = F ab F cd g ac g bd = F ab F cd (q ac - n a n c ) (q bd - n b n d ) = F ab F cd q ac q bd - 2F ab F cd n a n c q b 

2 



acbd 



IP IP „ ac r,' 
ra b r c dq q 



(A b - d b (A a t a ) - N a F ab ) (A d - d d (A a t a ) ~ N c F cd S j 



hd 



(4) 
(5) 



where A b — CtA b — t a d a A b + A a d b t a , and the action takes the form 
Sm = A A x^F ab F ab = -JL J dt d 3 xN^F ab F ab 



' f dt [ d 3 xN^/q~ 2 



16?r 



dt / d d 



s, 



TV 2 



(A b - d b (A a t a ) - N a F ab ) (Aa - d d (A a t a ) - N c F cd ) q bd + F ab F cd q ac q bd 



8ttN 



A b - dt (A a t a ) - N a F ab \(A d -d d (A a t a ) - N c F cd q 



16tt 



F ab F cd q ac q 



hd 



(6) 



r 



It follows that the conjugate momentum 7r a to the con- 
figuration variable A a is given by 



SS M _ Vq 



5A e 4tt7V 



(A d - 8 d (A a t a ) - N c F cd ) q ed , (7) 



measured by an observer with 4- velocity n a . Now the ac- 
tion can be expressed in terms of the canonical variables 
A„ and 7r a , 



which is a densitized vector field because of the presence 
of y/q. Its physical interpretation is as the electric field 



Sm {A 



a ,^) = jdt^ d 3 x (^L^ b q ab - ^F ab F cd q^q bd ^ . (8) 

by writing the integrand 



We can cast the action in equation (jSJ) into the desired form Sm = f dt J* s d 3 xn a A a — Hm 
in the following manner: 



Sm (A a , tt") 



dt I d 3 x 
= / dt I d 3 x 
= I dt I d 3 x 



h ( 2tt b yfq bd 

— ^7T tt q ab - I\ — 7T 7T q ab + — — b ah b cd q q 
yfq \y/q 16tt 

7T a (A a - d a (A d t d ) - N c F ca }j - N (^Ln a n b q ab + ^F ab F cd q ac q bd 
ir a A a + (A d t d ) d a ir a - N c ir a F ca - N (^n a n b q ab + ^-F ab F cd q ac q bd ^j 



(9) 



having integrated by parts in the second term. This com- equations of motion from equation @. First, since the 
pletes the Legendre transform and we can read off the 
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momentum conjugate to the time component of A a is 
absent, extremization of the action with respect to A a t a 
results in 



G = d a 7: a = 



(10) 



as the usual Gauss constraint. The total Hamiltonian of 
the Maxwell field then is 



H M = I d 3 x 



(A d t d ) d a n a + NVF a 



V I ^ n ~a~b n i jp jp „acJad 

N — — 7r ir q ab + ——b ab b cd q q 
\^q 16?r 



r 



(ii) 



with two contributions 



D c = n a F ca (12) 



and 



H = ^n a 7T b q ab + y^-F ab F cd q ac q bd (13) 

y/q 107T 

which, when added to the gravitational Hamiltonian, give 
matter contributions to the diffeomorphism and Hamil- 
tonian constraint, respectively. From (|13p we obtain the 
usual expression J d 3 xH for the energy of the electro- 
magnetic field. 



B. Equation of state 

The evolution equations can be obtained by evaluat- 
ing the Poisson brackets of A a and 7r a with the Hamil- 
tonian. Although we will not need the explicit form of 



J 



these equations, we present them in Appendix [Al for the 
sake of completeness. Here we rather determine energy 
and pressure from our canonical expressions (see also Q) 
in order to formulate the equation of state. The matter 
Hamiltonian is directly related to energy density [28| by 



1 SH 



A I 



V? SN 

and thus, from equation (|lip. it is 



2-7T 1 

^a^b„ i TP TP „ ac „bd 

P = — 7T 7r q ab + -—b ab b cd q q . 
q Voir 



The canonical formula for pressure is given by 



P 



SH 



-.q a b- 



AI 



.SH 



A I 



3N^/q wu Sq ab 3NJq H Sq ab 
as shown in Appendix [Bl This gives 



(14) 



(15) 



(16) 



P 



2 ,f } (^ a n b (q ab qef ~ 2q ae q bf ) + ^q™F ae F cf - ^F ab F ab q ef 



3AV9 VV? 



32tt 



( 



nN 



„ajb„ , V<jN b 

n 7T q ab + — — b ab b 
6Zn 



3N^q\^q 



2tt 1 

—a„b„ i ip rpab 
7T 7T q ab + ——b ab b 

q LOTT 



r 



(17) 



Finally, the equation of state can easily be obtained from 
([TBI and El): 



P 1 
w= — = - 
P 3 



which is the standard result. 



(18) 



III. QUANTIZATION 

Being interested in effects from quantum gravity, we 
have to quantize metric components in the matter Hamil- 



tonian (fTTj) , not just the matter field itself. Metric factors 
are thus not treated as a given classical background but 
become operators. This requires a quantum representa- 
tion, which can only be found if we also use canonical 
variables for the geometry. We thus need to use mo- 
menta of q ab even though they do not appear in the mat- 
ter Hamiltonian. 

In loop quantum gravity, the basic objects appropriate 
for a canonical quantization are constructed from a densi- 
tized triad Ef and the SU(2)-connection A\ = T l a - jK l a 
where T a is the spin connection compatible with the 
triad, K a the extrinsic curvature and 7 is the Barbero- 
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Immirzi parameter [j| [13, [Tl[ . Instead of the spatial met- 
ric q a b we thus use the densitized vector fields Ef which 
are related to the metric by EfE\ — q ab det q c d. 

These fields cannot be quantized directly but must be 
integrated suitably to remove local divergences in delta 
functions. The basic ingredient of a loop quantization is 
to use holonomies h e (A) = "Pexp J A l a e a Tidt £ SU(2) for 
all curves e c S and fluxes Fs(E) = J s Efn a T l d 2 y for all 
surfaces ScS where Tj are Pauli matrices, e a is the tan- 
gent vector to the edge e and n a the co-normal to the sur- 
face S. Thus the canonical quantization is performed by 
using holonomies and fluxes as operators, turning their 
Poisson brackets into commutators [12, [Tj] ■ A quantum 
representation is easily constructed by using states which 
are functionals of connections. Since holonomies are our 
basic connection dependent operators, they serve to gen- 
erate all states from a basic one which is just a constant 
on the space of connections. All states are then function- 
als depending on connections through holonomies, and 
they can be associated with graphs collecting the edges 
of holonomies used in the generation process. An or- 
thonormal basis can be determined explicitly in terms of 
spin network states [l4[ . 

An immediate consequence of this quantization is that 
fluxes and s patial g eometrical operators such as area and 
volume [lH fill [l7j have discrete spectra containing zero. 
Hence, their inverses do not exist as densely defined op- 
erators. However, a quantization of the matter Hamilto- 
nian such as demands the quantization of such in- 
verse expressions since, e.g., or the metric q a b which 
can only be obtained by inverting the densitized triad, 
appear in the matter Hamiltonian. Therefore, the quan- 
tization of the matter Hamiltonian seems, at first, to be 
seriously problematic. However, a well defined quantiza- 
tion is possible after noticing that the Poisson bracket of 
the volume with connection components, 

|4,y V\d^E\d 3 x} = 2v 1 Ge^ k e abc ^L = ^Ge l a 



ambiguities. However, several characteristic effects occur 
for any quantization choice such that they can be studied 
reliably with phenomenological applications in mind. 



V^det^ 

(19) 

amounts to an inverse of densitized triad components 
[lH . This is written here in terms of the co-triad e* 
from which we can directly obtain the metric q a b = e* e£. 
Similar expressions allow one to include the inverse de- 
terminant of the metric as we need it in the Maxwell 
Hamiltonian. The left hand side of p^|) does not re- 
fer to inverse densitized triad components and can be 
quantized: we can express the connection component 
through holonomies, use the volume operator and turn 
the Poisson bracket into a commutator. This observation 
enables us to quantize inverse powers of the densitized 
triad. Leading to well-defined operators, this quantiza- 
tion process implies characteristic modifications of the 
classical expressions such as (|TT1) on small scales, where 
densitized triad components are small. Moreover, since 
there are many different but classically equivalent ways 
to rewrite expressions like (1191) for which the quantiza- 
tion would give different results, there are quantization 



A. Effective Maxwell Hamiltonian 

Hamiltonian operators of a quantum theory can, in 
semiclassical regimes, be approximated by effective ex- 
pressions which amend the classical ones by quantum cor- 
rection terms. The general procedure, detailed in (l9ll20j. 
requires one to evaluate expectation values of the Hamil- 
tonian in suitable semiclassical states. A crucial ingredi- 
ent in loop quantum gravity is the discrete, non-local na- 
ture of states written in terms of holonomies as basic ob- 
jects. Although Hamiltonian operators on such discrete 
lattice states are quite complicated, expectation values 
can often be evaluated explicitly in perturbative regimes 
where one assumes the geometry to be close to a sym- 
metric one. This is certainly allowed in our applications 
to derive the effective equation of state of radiation in a 
flat FRW universe. The background symmetry implies 
the existence of three approximate spatial Killing vector 
fields Xj generating transitive isometries. We will only 
make use of this translational symmetry, not of the ad- 
ditional rotations in the construction of states. These 
vector fields can be used as a tangent space basis, thus 
denoting tensor indices for components in this basis by 
capital letters I,J,... 

The background symmetry also has implications for 
the selection of states of the quantum theory. A general 
quantization has to consider arbitrary states, but for ef- 
fective equations one computes expectation values only 
in states suitable for a semiclassical regime. For pertur- 
bative inhomogeneities, one can restrict lattices as they 
occur in general graphs to regular cubic ones and thus 
simplify geometrical operators. This has been developed 
; recently in 21] for metric perturbations as well as for 
a scalar field, and we can directly apply the same tech- 
niques to the Maxwell Hamiltonian. We refer the reader 
to this paper for more details. 



B. Gravitational variables and lattice states 

In a perturbative regime around a spatially flat 
isotropic solution, one can choose the canonical vari- 
ables to be given by functions (J) (x), kj(x)) which de- 
termine a densitized triad by Ef = p( % \x)5f and ex- 
trinsic curvature by K l a = k(i\{x)8 % a . Thus, one can di- 
agonalize the canonical variables compared to the gen- 
eral situation where all matrix elements of Ef and K l a 
would be independent. As seen in many symmetric mod- 
els, this simplifies the calculations considerably: it allows 
one to replace involved SU(2) calculations by much sim- 
pler U(l) calculations [H, [23[. SU(2) matrices arise be- 
cause loop quantum gravity is based on holonomies h e = 
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V exp(J e dte a A^Ti) of a connection A l a related to extrin- 
sic curvature. For unrestricted connections, holonomics 
can take any SU(2) value, but a diagonalization implies 
that all quantities can be reduced to a maximal Abelian 
subgroup U(l). Matrix elements of Hamiltonians and 
other operators can then be computed in explicit form. 

Using properties of the general loop representation 
mentioned before, basic variables of the quantum the- 
ory are, for a chosen lattice, U(l) elements i] v j attached 
to a lattice link e v j starting at a vertex v and pointing 
in direction Xf, and their conjugate fluxes F v j. The 
U(l) elements r/ v j appear as matrix elements in SU(2) 
holonomies h v j — Rer) v> i + 2rilmr] v j along edges e v j. 
Following the construction of the Hilbcrt space using 
holonomies as "creation" operators by acting on a state 
which is constant on the space of connections, a general 
state is a functional | . . . , (a v ,i, ■ ■ •) = Yl v / Vyi* • Allow- 
ing all possible values of assignments of integers fi v j € Z 
to the lattice edges e v j, this defines an orthonormal ba- 
sis of the Hilbert space. Basic operators are represented 
as holonomies 

Vvj\ ■ ■ ■ , th>',J, ■ • •) = I • • • , Mv,J + !,•■■) (20) 

for each pair (v, I) where all labels other than fj, v j remain 
unchanged, and fluxes 

Fv,i\ ■ ■ ■ ,IJ>v>,J, ■■■) = 2irj£p(iJ, v> i + (Jtv,-i)\ ■ ■ ■ ,Hv',j, ■ ■ ■) 

(21) 

where £p = \JhG is the Planck length and a subscript —I 
means that the edge preceding the vertex v in the chosen 
orientation is taken. These and the following construc- 
tions are explained in more detail in [211 ]. 

Effective equations are obtained by taking expectation 
values of the Hamiltonian operator and computing a con- 
tinuum approximation of the result (similar to a deriva- 
tive expansion in low energy effective actions). The result 
is a local field theory which includes quantum corrections. 
This is done by relating holonomies 

%iJ = exp(i / dtjki/2) mexp(M 'yki(v+I/2)/2) (22) 

e„,i 

to continuum fields fc/ through mid-point evaluation on 
the edges e v j (denoted by an argument v + 1/2 of the 
fields), and similarly for fluxes 

F v ,i= f p I (y)d 2 y^£ 2 p I (v + I/2). (23) 
Js Vtl 

Although the non-local basic objects do not allow us to 
define continuum fields at all spatial points, in a slowly- 
varying field approximation the mid-point evaluations are 
sufficient to define the continuum fields by interpolation. 
Here, £o is the coordinate length of lattice links. It does 
not appear in the quantum theory which only refers to 
states and their labels fi v j. This is independent of coor- 
dinates and only makes use of an abstract, labelled graph. 
The parameter £$ only enters in the continuum approx- 
imation since it is classical fields which are integrated 



and related to holonomies and fluxes. These continuum 
fields, or tensor components p 1 and kj, must depend on 
which coordinates are chosen to represent them. For the 
situation given here, the combinations p 1 := t^p 1 and 
ki := ioki, as they appear in holonomies and fluxes eval- 
uated for slowly-varying fields, are coordinate indepen- 
dent. 

A further operator we can immediately define is the 
volume operator. Using the classical expression V = 

J d 3 x y /\p 1 p 2 p 3 \ « J2 V ?W\p 1 p 2 p 3 \ = J2 V \ZIpV VI we 

introduce the volume operator V — 0/=i \J \^v,i I 
which, using ([2T]) . has eigenvalues 

3 

= ( 27r 74) 3/2 E II ( 24 ) 

v 1=1 

This operator is not only interesting for geometrical pur- 
poses, but also for making use of the identity (fT!))) or, 
more generally, 

{A i a ,V;} = 4n 1 GrV v r - 1 ei (25) 

which gives inverse powers of the densitized triad for 
any < r < 2 often appearing in matter Hamiltoni- 
ans. When quantizing this expression using holonomies, 
the volume operator and a commutator for the Poisson 
bracket, we obtain 

For symmetry, we use both edges e v j and e„,_j touching 

~ (r) 

the vertex v along direction Xf. The operator B y v I is 
obtained by taking the trace in (|26[) and using h v j — 
Re??jj,j + 2T/Im77 Uj /, 

with 

c v.i = r (Vv,I + Vvj) and Svj = TriVvj ~ Vv,i) ■ 

2 At ' 

Such expressions can be used for the electric field part 
of (|lip where the metric factor to be quantized is 

lab _ e a e b ^O e a e l 

in terms of the volume V v w £f)\/q(v) of a lattice site. 
This can then be quantized, using ([2^)1 with r = 1/2, to 

j^yz = (ioV v ~ 1/2 e])(£ Q V- 1/2 ey) = C^C^ . (28) 
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Noticing that the momentum 7r a of the electromagnetic 
field is quantized, just as the densitized triad, by a flux 
operator H v j := J s d 2 yn a ir a w ^q7t / (w), the whole 
electric field term can be written as 

= 2w [ d 3 xN(x)^w a w b « 2tt ^ JV(u)^ M^ 6 
= 2ir Y N{v)-^-R v jIi v j 



which is then quantized to 

= 27r^iV(u)C'^ /2) C'^ 2) n,, J n„ > . 



(29) 



1 



<7/j 



B v jB v 



with the magnetic flux H^/ := J* g d 2 yn a B a w EqB 1 ^). 

Magnetic flux components B 0j j are quantized using U(l) 
holonomies of the electromagnetic vector potential along 
closed loops transversal to the direction 7: 



B 



ajaxe Av,ajJ,a K K 



j,k <Tj,CT K e{±i} 



For the magnetic field term in (fTTj) , at first sight, a 
different metric expression arises: y/qq ac q bd which also 
involves inverse components when expressed in terms of 
the densitized triad. The term appears different from the 
electric field term and could thus be quantized differently. 
However, noting 



achd 



F ab F cd q ac q 



B e B f 



CeabCfcdq <Z 



t eah B*Blq jd e ahd q- x 



2q- 1 q ab B a B b 



in terms of the magnetic field B a — e abc Ff, c shows that 
the metric dependence is the same as in the electric part. 
We thus expect the same metric operator and corre- 
spondingly the same quantum gravity corrections in both 
terms, although different ones are mathematically pos- 
sible owing to quantization ambiguities. The magnetic 
contribution to the Maxwell Hamiltonian then is 

H B = — [ d 3 xN(x)^B a B b ^ —Y N(v)£^B a B 
8tt J y/q 8lr „ V9 

I 



We use the symbol A to distinguish an electromagnetic 
holonomy A from a gravitational one, r\. The loop holon- 
omy X v ,±j,±k is then computed around an elementary 
lattice loop starting in v in direction ±Xj and returning 
to v along iX^. Summing over J, K and the two sign 
factors a j and uk accounts for all four loops starting in 
v transversally to e v j. The resulting quantized magnetic 
part of the Hamiltonian is 



-E 



N(v)C { ^ 2) C^j 2) B v jB v j (30) 



with the same gravitational operator C^j^C^j as in 
the electric term. It is thus natural to use the same quan- 
tum operators and corresponding corrections in both 
terms, even though mathematically it is possible to quan- 
tize them differently. This aspect will be used in the 
following calculations. 



C. Effective Hamiltonian and equation of state 



As in |2l| we can include effects of the quantization of metric coefficients by inserting correction functions in the 
classical Hamiltonian which follow, e.g., from the eigenvalues [2l[ 

Ct /2) ({/V,J'}) = 2(27T 7 4r 1 / 4 | A V, / + ^-jfl^K + /V-id V4 (Kx + l*v,-K + 1| 1/4 - \l*v,K + l*v,-K ~ 1| 1/4 ) 

(31) 

(where indices J and K are defined such that ej jk 0) of operators C^j . Although for large fi v j these eigenvalues 
approach the function 



\/J>v,I + Hv,-l\\Hv,J + [l>v,-J\ 



expected classically for qij/^/q = \J\p x P 2 P z \lp l ' p J with 
a densitized triad Ef = p^Sf and using the relation (|2"Tj) 
between labels and flux components, they differ for values 
of Hvj closer to one. This deviation can, for an isotropic 



background, be captured in a single correction function 

a v , K = \y / ci i p({» v ,, r }f. 



v,I \ fJ>v,—I ) 



rij=i \Z\vv,j + vv,-j\ 

which would equal one in the absence of quantum correc- 
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tions. This is indeed approached in the limit where all 
AVj ^ 1) but for any finite values there are corrections. 
If all fj, v j > 1 one can directly check that corrections are 
positive, i.e. a v _K > 1 in this regime. Expressing the la- 
bels in terms of the densitized triad through fluxes (|2"Tj) 
results in functionals 



a = /3, the equation of state w can easily be computed 
and is modified as 



1 dloga 

w eS = - - 2— 

3 d log q 



(38) 



(33) 



which enter effective Hamiltonians. The general expres- 
sion one can expect is thus 



This modification is independent of the specific matter 
dynamics as in the classical case, and it results in an 
equation of state which is linear in p, but depends on the 
geometrical scales (and the Planck length) through a. 



H cS = / d xN 



2tt 

a[q c d]—n a TT b q ab 



0[q cd }^-F ab F cd q ac q bd 

l07T 

(34) 



with two possibly different correction functions a and 
depending on the lattice values p v ,i- As shown before, 
the case a = f3 is preferred, and we will see soon that this 
has implications for the effective equation of state. (In 
[J a Hamiltonian operator was introduced which did not 
use the same quantizations for metric coefficients in the 
electric and magnetic parts, thus giving a ^ (3. A quan- 
tization as described here, using the same quantization 
in both parts, was formulated in [24|. Phenomenological 
implications of a quantization of the latte r ty pe, con- 
cerning Lorentz invariance, are discussed in [25l| . ) There 
are other possible sources for corrections, such as higher 
order powers and higher derivatives of the electric and 
magnetic fields. But these terms would not be metric 
dependent and are thus not crucial for the following ar- 
guments. 

Now using (|34[) . we get the modified expression 



IV. COSMOLOGICAL APPLICATIONS 

In an isotropic and homogeneous universe (FRW), it 
follows from the FRW metric and Einstein's equation 
that the evolution of the energy density is given by the 
continuity equation, i.e., 



p + 3-(p + P) 
a 



0, 



(39) 



where a is the scale factor and the dot indicates a proper 
time derivative. Using the definition of the equation of 
state and eliminating the time derivative, this equation 
can be cast into the following useful form: 



dk)gp(a) 
d log a 



-3(l + w(a)) 



(40) 



Here we have shown the dependence of the equation of 
state on the scale factor explicitly. It can easily be shown 
that the solution to the above equation is 



1 ab 6H M q a bSH M it d , , 6 , 

= —pn 7r q cd [a + 2q da/dq ) 



N q Sq ab 



p(a) = p exp 



N 5q ab ^fq 
+ ^ F cdF Cd (f3 + 2q ab S/3/Sq ab ) 



(1 + w(a)) dloga 



(41) 



(35) 



depending on a and /3. For a nearly isotropic back- 
ground geometry, for instance, a only depends on the 
determinant g of the spatial metric and, from Appendix 
B, q ah 8a/5q ab = — 3gda/dg, which we assume in what 
follows. 

The modified energy density and pressure then are 

Pes = —K a K h qabU + ^-F ab F cd q ac q M p (36) 
q l07T 

2tt 

3P cS = — TT a 7r b q ab (a - 6gda/dg) 

q 

+ -^F ab F ab 09 - 6gd/3/dg) 

1D7T 

271" a b f, r AlQga 

= — 7r 7r q ab a 1 - 6— 

g V dlogg 

+ ±F*F*0 (l-Bp£). (87. 
167T \ dlogg ' 



where po is the integration constant. Now by inserting 
the modified equation of state in the radiation era, (|38|) 
with q = a 6 , we obtain 



p(a) = p a(a)a~ 



(42) 



It follows easily from ([33)) , (|3"6")l and ([3T)l that the classical 
behavior is reproduced for a = (3 = 1. Interestingly, for 



Again, for a = 1, we retrieve the classical result p(a) oc 
a -4 . Therefore, loop quantum gravity corrections in- 
duced by discreteness of the flux operator are reflected 
even in the evolution of the FRW universe. 

Although one can write a as a function of the scale 
factor for perturbations around a flat isotropic model, 
it is important to note that corrections are well-defined 
even though the scale factor can be rescaled arbitrar- 
ily. One can express a as a function of a only after 
coordinates have been specified, such that there is no 
ambiguity in relating the scale appearing in a (such as 
the Planck length) to a. More precisely, we first ob- 
tained a[p / (w)] in (|3"3"|) with lattice values for the fluxes 
p 1 (v) = iyp 1 (v) which are coordinate independent while 
p 1 would be rescaled just as the scale factor. The quan- 
tum state, through its lattice building blocks, unambigu- 
ously determines the magnitude of the elementary vari- 
ables as they appear in corrections. Under rescalings or 
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other coordinate changes, both the classical field p 1 and 
the coordinate form of the lattice change in such a way 
that elementary fluxes remain unchanged. In a nearly 
isotropic context, for instance, one has p 1 (v) sa p = i^a 2 
spatially constant which can be related to the Hubble 
scale by Np 3 / 2 = H~ 3 . Here, we use the number N of lat- 
tice sites of elementary area p in the Hubble volume H ~ 3 
as a measure of how fine the lattice is. Inserting all this 
in correction functions yields a(p) = a(N~ 2 / 3 H~ 2 ) ex- 
pressed purely in terms of coordinate independent quan- 
tities. The function N enters as an additional ingredient 
to describe the microstructure of the underlying quantum 
state. In a given solution including the time dependence 
H(t) of the Hubble scale as well as a function N(t) de- 
scribing the quantum state one could relate all this, in 
a secondary step, to the scale factor a(t). But since the 
scale factor is not the primary argument of correction 
functions, there is no problem with rescalings. See also 
(26| for further clarifications of this issue which was not 
clear in all the literature on purely homogeneous models. 



V. DISCUSSIONS 

We have derived here the equation of state of the 
Maxwell field in a canonical form, including corrections 
expected from loop quantum gravity. In the canonical 
derivation, the reason for a linear equation of state, which 
is trace-freedom in the Lagrangean derivation, is the fact 
that the same metric dependent factor q a b / ^/q multiplies 
both terms in the Hamiltonian. The Maxwell Hamilto- 
nian is thus simply rescaled if the metric is conformally 
transformed, which explains the conformal invariance of 
Maxwell's equations. This is special for the Maxwell 
field and different from, e.g., a scalar field with a non- 
vanishing potential. 

The same fact allows one to quantize the Hamiltonian 
in a way which affects both the electric and magnetic 
term in the same way, at least as far as the metric de- 
pendence is concerned. One then obtains a single cor- 
rection function a = j3 which only corrects the metric 
dependence of the total scale of the Hamiltonian. In this 
sense, conformal invariance is preserved even after quan- 
tization. (But this would not be the case if a quantization 
is used which results in a ^ (3.) 

This preservation of the form of the Hamiltonian ex- 
plains why we are still able to derive an equation of state 
independently of the specific field dynamics and that it 
remains linear. However, the classical value w = ^ is 



corrected due to quantum effects in the space-time struc- 
ture. This modification is also understandable from a 
Lagrangean perspective, together with basic information 
from the loop quantization. Employing trace freedom of 
the stress-energy tensor to derive the equation of state, 
we have to use the inverse metric in g ah T a b- But from 
loop quantum gravity we know that, when quantized, 
not all components of the inverse metric agree with in- 
verse operators of the quantization. For the scale factor 
of an isotropic metric, for instance, we have a -1 ^ "a -1 " 
since the right hand side is not even defined @. While 
the left hand side is defined through identities such as 
(|19p. it satisfies a~ l a ^ 1 and thus shows deviations from 
the classical expectation a~ l a = 1 on small scales which 
were captured here in correction functions. As derived in 
detail, this implies scale dependent modifications to the 
equation of state parameter w e f{. 

The result can also be interpreted in more physical 
terms. The classical behavior p{a) oc a -4 can be under- 
stood as a combination of a dilution factor a~ 3 and an ad- 
ditional rcdshift factor a~~ l for radiation in an expanding 
universe. As we have seen, this is corrected to a(a)a -4 
where a(a) corrects the metric factor q a b/^/Q ~ a 5 a },. 
Since this is only a single inverse power of a for an 
isotropic solution, we can interpret the result as saying 
that only redshift receives corrections due to quantum ef- 
fects on electromagnetic propagation. The dilution factor 
due to expansion is unmodified, except that the back- 
ground evolution a(t) itself receives corrections. This 
agrees with the result for dust, which is only diluted and 
has an unmodified equation of state even after quantiza- 
tion [29j. Unlike dust, for radiation one has to refer to 
the inhomogeneous field and its quantum Hamiltonian to 
derive a reliable equation of state, as presented here. 
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APPENDIX A: EQUATIONS OF MOTION 

It is straightforward to derive the equations of motion 
for the canonical variables A a and ir a from the Poisson 
brackets of each of these variables with the matter Hamil- 
tonian Hm- Then 



A a = {A a , H M } - ^ = d a (A c t c ) + N c F ca + ^n c q ca , 
ou a y/q 



and 



7T Q = K,ff M } 



r tt 

— ™ = 8 C (N c Tr a ) - d d (N a n d ) - 4d c (N^F ef q ec qf a ) 
0A a 



(Al) 



(A2) 
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The modified Hamiltonian gives rise to the following new set of equations of motion: 

A a = {A a , H eS } = = d a (A c t c ) + N c F ca + —a (q) 7r c q ca , 



(A3) 



and 



TT a = {TT a , H cS } 



' SA a 



8 C ( W) - 8 d (N a n d ) - 48 c (N0 (q) ^F ef q ec qf a ) , 



r 



(A4) 



where H e g is the effective Hamiltonian of the Maxwell's 
field (H M with a and fi inserted). 

APPENDIX B: PRESSURE 

The general, thcrmodynamical definition of pressure is 
the negative change of energy by volume, which we can 
write as 



P = - 



1 SH 
N5^ 



(Bl) 



whenever the Hamiltonian H = J d 3 xN(x)H(x) is de- 
pends isotropically on the metric. Otherwise, one has to 
use all components of the stress tensor SH/6q ab which is 
not proportional to the identity. The derivative by the 
determinant of the metric can be expressed in terms of 
metric components by using a suitable change of vari- 
ables which includes q as an independent one. We thus 
introduce q a b =: q^^qab with Aetq ao = 1 such that 
dq a b/dq = ^q~ 1 q a b where all components of q a b are kept 
fixed in the partial derivative. This is exactly what we 
need to compute pressure since only the volume but not 
the shape of the fluid is varied. This change of variables 
implies 



6 or 6 



and thus 



'■lab 



ab 



dq 6q ab 



ab 



P = — 



SH 



q a b- 



iN^/q 8q ab 



Sq a b 



(B2) 



We can also verify this by comparing the dynamical ef- 
fects of H on the metric with the Raychaudhuri equation 
expressed in terms of the canonical variables which for 
simplicity we do for homogeneous metrics. Using the 
following definitions for the extrinsic curvature tensor 
K a b = V a n{, (which turns out to be automatically spa- 
tial and symmetric without projection if homogeneity is 
used), the expansion parameter 9 — K a bq ab and the shear 



Oab 



(ab) 



8q a b, the canonical momentum conjugate 



to q a b derived from the gravitational Lagrangian is 



ab 



16ttG 



where G is the gravitational constant. Then the Ray- 
chaudhuri equation in terms of the canonical variables 
takes the following form: 



(B3) 



The canonical equations of motion, in the presence of a 
matter Hamiltonian H added to the gravitational Hamil- 
tonian to form -ff Total, become 



^Total 167TG7V 

lab = ,_ nh = 7=— (27r a6 - q ab ir c ) + 2D {a N b) 



5ir ab 



(B4) 



and 



.ab = J*** (S )R ab h 3)R ab\ + ^ ab ( 

Sq ab 16irG V 2 J y/q \ 



1 2 \ 32ttGN b i _„,. / !._ „,, 



d7T ~r )-^r q v nc ~^ 



t— + (D a D h N - q ab D c D c N) + y/qD c ( ) - 2^ a D c N b \ (B5) 

Sq ab 16irG V V? / 



I 

where D a is the derivative operator compatible with q a b- function N yields the Hamiltonian constraint equation 
Variation of the total action with respect to the lapse 
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Upon inserting equations (|B4|) , (|B5[) , and (|B6|) into equa- 
tion (|B3|) . the Raychaudhuri equation becomes 



Raychaudhuri equation is found to be 



1 - -Itfl-rr^rr ^ 

N 3 aab N^q" ' Ny/q^Sq* 



ff + ^^gah— (B7) 



D"D„ V S; : (V/;, [ ^1) + l ^^a DcNa 



which, for a homogeneous universe, reduces to 



-e 2 - a ab a ab - AwG (p + 3P) . 



(B9) 



N 



Now comparing equation ()B8|) with equation (|B9[) . we 
verify the canonical formula for the average pressure for 
a perfect fluid distribution in an anisotropic geometry 



l = - 1 -e>-*°»* ab -^H + ^ q ™ (B8) 
N 3 N^q N^q qab Sq ab V ' 



2 <57J 

= <7ab 



5H 



(BIO) 



On the other hand, for a perfect fluid distribution, the 
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